Generalized DMPK equation for strongly localized regime - numerical solution 
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Generalized Dorokhov-Mello-Pereyra-Kumar (GDMPK) equation [K. A. Muttalib and 
J. R. Klauder, Phys. Rev. Lett. 82, 4272 (1999)] has been proposed for the description of the 
electron transport in strongly localized systems. We develop an algorithm for the numerical solu- 
tion of this equation and confirm that GDMPK equation correctly describes the critical and localized 
regimes. Contrary to the original DMPK equation, the generalized one contains also an information 
about the dimension of the system. In particular, it distinguishes between the two and the three 
dimensional models with the same number of transmission channels. 

PACS numbers: 73.23.-b, 72.10. -d, 02.70.Ns, 05.10.Gg 



Two decades ago, there was the belief that the dis- 
tribution of the logarithm of the conductance, hig, of 
the strongly disordered electron system is Gaussian, in- 
dependently on the dimension of the system. This 
paradigm was based on the two-terminal expression for 
the conductance, [H 
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Eq. ((T|) expresses the conductance in terms of parameters 
Xi which determine the eigenvalues of the matrix vt (t is 
the N x N transmission matrix and ./V is the number of 
open channels). Since all parameters Xi should increase 
linearly with the system size, x% oc L, it is natural to 
expect that in the limit of large L, only the contribution 
of the first channel (that one with the smallest parameter 
x) survives and In 5 rs —x\ + In 4. In the analogy with 
the one-dimensional (ID) case (N — 1), the distribution 
p(lng) should be Gaussian. This is consistent with the 
solution of the Dorokhov-Mello-Pereyra-Kumar (DMPK) 
equation, [2J for very long quasi-one dimensional (Q1D) 
weakly disordered systems. [3j 

However, numerical data [4[ for strongly disordered 3D 
samples showed that p(ln g) is not Gaussian. Although 
all Xi oc L, the mean values of the differences Xi + ± — 
Xi, i = 1,2,... are constant, independent on the size 
of the system. Therefore, higher parameters, X2, Xa, ■ ■ ■ , 
affect the form of the distribution p(ln g) in 3D. Since the 
deviation from Gaussian distribution cannot be obtained 
from the DMPK equation, the generalization of the last 
has been proposed by Muttalib and co-workers. @, Q 

The derivation the DMPK equation^ is based on the 
parametrization of the transfer matrix T. For spin-less 
electrons and the time reversal symmetry of the system, 
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where A is a diagonal matrix with diagonal elements 
A, = (1 + cosh2xi)/2, and u, v are unitary matrices. 



0] For sufficiently long systems, it is assumed that the 
elements of matrices u, v and A are statistically inde- 
pendent. Contrary to the classical DMPK, the general- 
ized DMPK equation (GDMPK) contains additional 
N(N — 1) parameters, 



Ki 
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where (...) means a mean value over the realization of 
the disorder. The GDMPK equation reads 
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Here, P({x}) is the common probability distribution of 
all parameters x, £ is the mean free path, I. The inter- 
action potential f2({x}) reads 



£l({x n }) = - J2i <3 Iji m I sinh (xj) - sinh 2 (xi)| 



X;il 1 ln|sinh(2x i )|, 
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and 7y = 2K lj /K ii , 

Parameters Kij depend on the strength of the disorder. 
For weak disorder, — [1+%]/ (N+l), 7y = 1 and the 
GDMPK equation reduces to the classical DMPK equa- 
tion. [2j, |8| For general disorder we expect that K con- 
tains an information about the strength of the disorder, 
the top olog y and the dimensionality of the disordered 
system [lQj . For instance, d Ku ~ 1 and K 12 - l/L 
in the localized 3D regime. We expect that the transi- 
tion from the metallic the insulating regime is associated 
with (continuous) change of parameters K [^, [l(| . and 
the GDMPK equation with correct choice of parameters 
K describes the transport in the metallic, insulating and 
even in the critical regime. Note, the classical DMPK 
equation is not applicable to strongly disordered systems 
and contains no information about the dimension. 

In Refs. 0, [T3|, a simplified GDMPK equation was 
introduced, with only two parameters: the diagonal ele- 
ments Ku = K\\ for all i, and the off-diagonal elements 



2 



Kij = K\2 for all i ^ j. [ll| An approximate solutions 
of such GDMPK equation, based on the saddle-point 
method [13] was obtained in Refs. d, and the solu- 
tion was compared with numerical data obtained by the 
transfer matrix method (TM). Although quite satisfac- 
tory, agreement between approximate solutions and nu- 
merical data was obtained, a quantitative analysis of the 
solution of the GDMPK equation is still missing. Con- 
trary to the DMPK equation, which is exactly solvable 
[HI , no exact solution of the GDMPK equation is known. 
It is therefore highly desirable to solve the GDMPK equa- 
tion numerically. 

In this paper, we present an algorithm for the numer- 
ical solution of the GDMPK equation. Our method is 
based on the mapping of the GDMPK equation onto the 
Langevin equation, which describes the diffusion of N 
particles interacting with potential Q({x}), given by Eq. 
(|5|) . Simulating such diffusion, we obtain numerical solu- 
tion of classical DMPK equation and the GDMPK equa- 
tion. These solutions are compared with the numerical 
TM data for the tight-binding Anderson model, 
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which describes the transport of a single electron on the 
d-dimensional lattice. Random energies e r have zero 
mean value and variance 1/12, W measures the strength 
of the disorder. Hopping term between two nearest neigh- 
boring sites r and r' is tn = 1 in the direction of the 
propagation and t± = 0.4 in the perpendicular direc- 
tions. For this choice of parameters, there are no evanes- 
cent (closed) channels in 2D and 3D systems. 14] In 3D, 



model ([5]) exhibits the disorder induced transition from 
the metal to the insulator at the critical point W c ~ 9.2. 

Our numerical method uses the fact that the GDMPK 
equation Q is a special Fokker-Planck diffusion equation 
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for the one dimensional diffusion of N particles located 
at x in "time" s = L/t. The common probability distri- 
bution P({x}, s) determines the positions of all particles 
at time a. In Eq. ((7J, D^' and D^' are the drift and 
the diffusion coefficients. The Fokker Planck equation 
([7]) describes the random process given by the Langevin 
equation, [lj| 



dxi/ds = hi({x},s) +# i ({a:})r i (s), 



(8) 



with random white noise force T(s), (r(s)) = and 
(r(s)r(s')) = 2<S(s - a'), and coefficients 
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FIG. 1: The test of the numerical algorithm. The classical 
DMPK equation is solved for N = 9 channels. Results are 
compared with the TM data for the Anderson model on the 
lattice 3 x 3 x L and disorder W = 2. Solid lines are the 
TM data for L — 30 (used as initial condition for the DMPK 
equation), and L = 50, 70, and 100. Symbols show solutions 
of the DMPK equation for s — 1.8 (circles), 3.6 (triangles) 
and 6.35 (squares). For these parameters, both the TM and 
the DMPK solutions have the same value of the mean con- 
ductance (g). 



Comparison of Eqs. ([!]) with Eq. ([7]) gives 
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We simulate the Langevin force T, integrate the equa- 
tion of motion Q and take the average for a large number 
of realizations. The numerical integration of the equation 
© gives 
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where N is the number of particles and 5s is the "time" 
step. Statistically independent Gaussian variables Wi n 



have zero mean and variance (wj n Wkr. 
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ter n integration steps we obtain the stochastic variables 
%in = %i(a) at "time" s = n x Ss. 

In simulations, we have to keep the ordering, xi < 
X2 < ■ ■ ■ < xn- Therefore, the position of the zth particle 
is restricted as £;_i n < Xi n < Xi+in-i- To keep this 
constrain, the "time" step, Ss must be very small. For 
xin the left boundary is 0. To avoid numerical overflows, 
we introduce and we use the "cutoff" C, xnu = C,Vn. 

We are not able to start the numerical algorithm with 
the ballistic initial values, lim^^o P = J|- <5(Aj — + ), 
used in the formulation of the DMPK equation. To avoid 
this difficulty, we use the TM method and calculate first 
the distributions p(xi) for the Anderson model for a short 
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FIG. 2: (Color online) The 3D Anderson model 7 x 7 x L 
with disorder W = 29 (strongly localized regime). Solid 
line shows the TM data forL = 6, circles are the TM data 
for lengths L = 8 and L — 10, which are compared with 
the GDMPK data for s = 10.5 and 20, respectively. In the 
GDMPK equation, we use either entire matrix Kij (triangles) 
or the two-parameter model (squares) with K\\ — 0.5225 and 
Ki2 = 0.0244. Estimated mean free path £ ~ 0.2 is consistent 
with [lbT |. The GDMPK data for p a ii(a;) are shown by dashed 
and dot-dashed lines for model with complete matrix K and 
the model with Ku and Km, respectively. 




5 10 5 10 




5 10 15 20 25 -20 -15 -10 -5 



FIG. 3: (Color online) The 2D Anderson model 49 x L and 
W = 26. Initial condition are given by the TM data for L — 6 
(solid line) . The TM data for L = 8 and L — 10 are compared 
with the GDMPK data for N = 49 channels and s = 9.5 and 
19 (t ~ 0.21). The symbols have the same meaning as in 
Fig. [2] We used Ku = 0.712 and K 12 = 0.0837 in the two 
parameter GDMPK equation. 



length Lq of the system. We use these distribution as 
initial values for the GDMPK equation. The TM data for 
L = Lq provides us also with the values of parameters K . 
The solution of the GDMPK equation is compared with 



the TM result for length L > L . In the TM method, 
the length is denned as a number of lattice sites along 
the propagation direction. For a given L, we find the 
length s in the GDMPK equation from the condition that 
the mean conductance (or the mean of logarithm of the 
conductance in strongly localized regime) found by the 
two methods coincides. We find in this way also the mean 
free path (in units of the lattice period of the Anderson 
model): £ = (L — L )/s. 

Our results are summarized in Figs. 1-5. Compared 
are four distributions: p(x±) andp^) for the two small- 
est parameters x, the distribution p a n(x) of all parame- 
ters x, p a n(x) = S(x — Xi)), and the distribution of 
the conductance p(g) or of the logarithm of the conduc- 
tance, p(ln g). 

In order test the numerical algorithm, we solve first 
the classical DMPK equation. Fig. [1] shows the data for 
the weakly disordered Q1D systems with N = 9. The 
length of the system varies between L = 30 (initial con- 
dition for the DMPK) and L = 100. The "time" step 
used in the solution of the DMPK equation 6s = 10~ 4 
and 63.500 iterations were performed. With the use of 
appropriate length s, we obtained perfect agreement be- 
tween the DMPK and TM data, and also estimate the 
mean free path, £ « 11, consistent with I ss 9.2 estimated 
for the 3D anisotropic systems [Toj ] . 

In Figs. Hand[3]we present the solution of the GDMPK 
equation for N = 49 channels in the strongly disordered 
regime (disorder W — 26). In Fig. [5] we simulate the 
transport through the 3D system 7 x 7 x L. The TM 
data for L = 6 provides us with the initial distributions 
for the GDMPK as well as with values of parameters K. 
Data for L = 8 and L = 10 are then compared with 
two solutions of the GDMPK: in the first simulations, 
we use the entire matrix Kij, and in the second simu- 
lation, we substitute all diagonal elements by Ku and 
all off-diagonal terms by Ki2- Fig. [3] summarizes the 
TM results for the transport through the 2D strongly 
disordered system 49 x L (W — 26) compared with the 
corresponding solution of the GDMPK equation. 

Both Figures show that the GDMPK data for p(xi), 
p(x2) and p(\ng) agree with results of the TM simula- 
tions. Since the TM data were obtained for systems of 
different dimension, we conclude that the GDMPK equa- 
tions distinguishes between the 2D and 3D disordered 
systems. The information about dimension is given in 
parameters Ku and Kyi. 

The difference between the 2D and 3D systems is more 
clearly visible in Figs. [4] and [5] which compares 3D TM 
data for the critical disorder W c — 9.2 and 2D data for 
W = 8.6. Although both system have approximately the 
same mean free path, the conductance distributions differ 
from each other when L increases. 

Obtained results confirm that already the simplified 
two-parameter GDMPK equation determines correctly 
the statistics of the conductance, both in the localized 
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FIG. 4: (Color online) The 3D Anderson model 7 x 7 x L 
with critical disorder W c = 9.2 Solid lines show the TM data 
for L = 6 (initial condition). The TM data for L = 10 are 
compared with the GDMPK data for N = 49 and length 
s = 3.25 and 6.375 {£ » 0.63). 




FIG. 5: (Color online) The 2D Anderson model 49x1 with 
disorder VP = 8.6 Solid lines: the TM data for L = 6 (Initial 
condition). Circles: the TM data for L — 8 and L = 10 
compared with the GDMPK data for s = 3.375 and s = 6.5, 
respectively (£ ps 0.615). Symbols have the same meaning as 
in Fig. 2. 



and critical regimes. However, the total distribution of 
all parameters x, p a \i{x), is reproduced only for rather 
small values of x. The possible reason for this discrep- 
ancy is that the TM data for the entire matrix K might 
not be sufficiently accurate for high indices i and j, so 
that we do not know true values of repulsion constants 
7ij for higher channels. In the two parameter model, 
Kn and K12, might overestimate the repulsion of higher 
channels (for instance 71, ~ i -1 / 2 (10]). This overestima- 
tion broads the distribution Pa3i{x). Interestingly, the 2D 



data for p a i\(x) are more accurate than the data for 3D. 
Nevertheless, as argued previously 17j . and confirmed 
by our data shown in Figs. 2-5, we do not need to know 
the distribution of all parameters x for the description of 
the statistics of the conductance. Only small portion of 
channels is important for the description of the critical 
regime in 3D. 

In conclusion, we presented numerical solution of the 
generalized DMPK equation. Our results confirm that 
the generalized DMPK equation describes correctly the 
electron transport in the localized and the critical regime. 
The information about the dimension of the system 
is g iven by the additional parameters, K\\ and K12 
1 11 ] . We hope that analytical solution of the general- 
ized DMPK equation could provide us with the analyt- 
ical description of the electron transport in the critical 
and localized regimes in various dimensions. 
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